< 






> 

On 

in 

o 



X 



Exact dynamics of one-qubit system in layered 

environment 



o 

(N 

m ; M.Mahdian a * ,H. Mehrabpour a 



department of Theoretical Physics and Astrophysics, University of Tabriz, Tabriz 51664, Iran. 

April 23, 2013 



"E-mail:Mahdian@tabrizu.ac.ir 



Abstract 

We investigate the exact evolution of the reduced dynamics of a one qubit system as 
central spin coupled to a femionic layered environment with unlimited number of layers. 
Also, we study the decoherence induced on central spin by analysis solution is obtained 
in the limit N — > oo of an infinite number of bath spins. Finally, the Nakajima-Zwanzig 
(NZ) and the time-convolutionless (TCL) projection operator techniques to second order 
are derived. 

Keywords : spin star model, decoherence, Nakajima-Zwanzing and time- 
convolutionless 

1 Introduction 

In quantum mechanics, quantum information is physical information that is held in the state 
of a quantum system [I]. Quantum information theory focuses on the amount of accessible 
information [2], it can be regarded as the theory for quantitative evaluation of the process of 
extracting information^ H]. Every quantum system encountered in the real world is an open 
quantum system and the theory of open quantum systems describes how a system of interest 
is influenced by the interaction with its environment. This interaction often leads to a loss of 
the quantum features of physical states and has a great impact on the dynamical behavior of 
the open system due to the non-unitary characteristic of the time evolution, although much 
care is taken experimentally to eliminate the unwanted influence of external interactions, there 
remains, if ever so slight, a coupling between the system of interest and the external world [3 
[6]. One kind of open quantum systems study to describe the information extraction process 
in quantum information are the spin star systems, where a central spin -| particle couples 
to a spin bath of N spin-| particles and they have attracted a vast amount attention in 
the quantum community [T]-[Tl] because they are of significance and of interest due to their 
high symmetry, strong non-Markovian behavior and also as one of the best candidates of the 



spin-qubit quantum computation[in]-[Tl]. This is even more relevant when environmental 
influences of a non-Markovian nature, such as those due to memory-keeping and feedback- 
inducing system-environment mechanisms, are considered [5]. 

The spin star configuration can also describe decoherence model [T5] because the coupling 
of an open quantum system with its environment causes correlations between the states of the 
system and the bath[5]. the correlations exchange the information between the open quantum 
system and its environment and the environment-induced, dynamic destruction of quantum 
coherence is called decoherence [16J. In the language of state and density matrix, the super- 
position of the open quantum systems states is destroyed after tracing over the environmental 
degrees of freedom and the systems reduced density matrix turns into a statistical mixture. 

Motivated by this consideration, in this paper, we consider layered environment with a spin 
at the center of layers to study a generalized spin star system which can be solved exactly. It 
must be noted that in the model, degeneracy for coupling coefficients are considered. 

The paper is organized as follows. In Sec. (2), we introduce the model investigated, a spin 
star model involving a Heisenberg XX coupling in Sec. (2.1), and determine the exact time 
evolution of the central spin in Sec. (2.2). Therefore, if we equalize all coupling coefficients 
with together, we obtain the result of Ref.pl)]. In Sec. (3) we assume two different layers 
of environment and compute our model with this assumption. Furthermore, we analyze the 
limit of an infinite number of bath spins, discuss the behavior of the von Neumann entropy 
of the central spin, and demonstrate that the model exhibits complete relaxation and partial 
decoherence. The non-Markovian approximation techniques are discussed in Sec. (4). In this 
section the dynamic equations found in the second order of the coupling are introduced. It is 
also demonstrated that the prominent Born-Markov approximation is not applicable to the spin 
star model. Of course, the Born-Markov approximation is second order Nakajima-Zwanzig. 



2 Exact Dynamics 

2.1 The Model 

We consider a spin star configuration which consists of N+l localized spin-| particles. One of 
the spins is located at the center of the star, while the others are on concentric circles with 
different radii surrounding the central spin, layer by layer the difference in radius is because 
the coupling coefficients between layers spins and the central spin are taken differently. It 
must be noted that in this model degeneracy for coupling coefficients is considered, because 
naturally some of spins are in relation with the central spin by a constant coupling coefficient 
which are located in one layer. By considering such model, the most general model of single- 
qubit spin-Star for fermionic particles with fixed fermionic environment is made, Fig(l). This 
model explains how particles of bath with different coupling coefficients can be used to control 
time of decoherence. Because degeneracy coefficient specifies the number of the particle layer, 
decoherence-time by each layer can be controlled with different degeneracy coefficient and 
we should not lose the effect of degenerate factor. However, we consider our model by the 
below description and explanation. The central spin a interacts with the bath spins a^ via a 
Heisenberg XX interaction [T7j represented through the Hamiltonian 

H = 2(a + ~_ + a_~ + ), (2.1) 

where E± are denoted as follows 



n 



E^ J +' (2-2) 



S_ = £;«„■£• (2.3) 

Here, a coefficients specify interaction between system and environment and is dependent from 



distance. Also we have 

N u 



J5N£°i (2-4) 




Figure 1: The figure depicts the general layered environment with one spin in the center of 
layers. 



here, (j, — 1,2, ... ,n for n different layers of bath. Also, we have 

that represents the raising and lowering operators of the jth bath spin. The Heisenberg XX 
coupling has been found to be an effective Hamiltonian for the interaction of some quantum 
dot systems |18j . Equation(l) describes a very simple time independent interaction with equal 
coupling strength a,\ for N\ of first bath spin and a 2 for N 2 of second bath spin to a n for N n of 
nth bath spin. It is invariant under rotation around the z—axis. The operator J = | X/«=i afl 
represents the total spin angular momentum of the bath (units are chosen such that H — 1). 
Therefore the central spin thus couples to the collective bath angular momentum. 

We introduce an Orthonormal basis in the bath Hilbert space Hb consisting of states 
\jn, m M , x) where /i is 1 to n. These states are defined as eigenstates of J3 (eigenvalue m) and 
of J 2 (eigenvalue j(j + l)). The index x labels the different eigenstates in the eigenspace Kj >m 
belonging to a given pair (j,m) of quantum numbers. As usual, j M < -*■ and — j M < m < j^ 
where n is 1 to n. The dimension of Kj >m is given by the expression [191 EQ] 
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We assume that the initial state of the composite system be a product state. That is 



p(0) = Ps(0)®p B (0). 



(2.6) 



We can calculate the reduced density matrix of the quantum system in the following expression. 



Ps(t) = tr B (Up(0)rf). 



(2.7) 



The above equation is obtained by doing partial-trace on bath and also U is the unitary 
operator which is defined as follows 

U = exp(-iHt). 

The reduced density matrix is completely determined in terms of the Bloch vector 



x(t) 



( wi(t) ^ 
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Tr(ap s (t)), 



(21 



through the relationship 



Ps(t) 



(2.9) 



1 + u 3 (t) LUi(t) - iu 2 (t) 
idi(t) - iuj 2 (t) l-u 3 (t) 

We note that the length q(t) = \x{t)\ of the Bloch vector is equal to 1 iff ps(t) describes a pure 
state, and the von Neumann entropy S of the central spin can be expressed as a function of 
the length q(t) of the Bloch vector: 



S = Tr[-p s \np s ] = In 2 - ^(1 - q) + ^(1 + q) ln(l + q). 



(2.10) 



The initial state of the reduced system at t=0 is taken to be an arbitrary (possibly mixed) 
state 
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(2.11) 



while the spin bath is assumed to be in an unpolarized infinite temperature state: 



p B (0)=2~ N I B . 



(2.12) 



Here, I B denotes the unit matrix in H B and N is N\ + N 2 + • • ■ + N n , and we have defined the 
u)± as linear combinations of the components u)i t 2 of the Bloch vector 



w± 



UJ\ ± IU2 



(2.13) 



2.2 Reduced System Dynamics 

In this section, we will derive the exact dynamics of the reduced density matrix ps(t) for our 
given model. We obtain the evolution of central spin with n different coupling coefficients that 
should be used from Eq.(7) until the solution model is exact. This yeilds 



PS 



(t) = **{<'£ E 7f ho 
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It can easily be verified that 



(2.14) 
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(2.15) 



and 



H 2k+i = 2 . 4 fcr s _ , a _ E _ 



(2.16) 



We note that such simple expressions are obtained since a term cr 3 J 3 is missing in the interaction 
Hamiltonian. We substitute the last two equations into Lindblad equation [in] as follows 



L l p = i l J2{-l) n 
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(2.17) 



to get the formulas 



tr B { L 2l+1 Ps (0) ® 2-^/^} = 



(2.18) 



and 



tr B { L 2l p s (0) ® 2""/*} = £(-4) J ^ ' ' 
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which hold for all 1=1,2,. . . . Here,we have introduced the bath correlation functions 

Mi = ^tr B {(E + E.) 1 }, 
Tt n = iir B {(S + S_)'-"(S_S + )"}, 



n ~ 2 N 

where we have product ^±E> T as follows 



(2.19) 

(2.20) 
(2.21) 






(2.22) 



Of course, Eq.(18) is zero because it is 



(j,m\J±\j,m) = 0. 

We will come back to these correlation functions when we discuss approximation techniques 
in Sec. (4). 

Using the formulas (18) and (19) in Eq.(14) we can express the components of the Bloch 
vector as follows, 

w±(«) = /±(<)w±(0), (2.23) 



u 3 (t) = f z (t)u 3 (0), 



(2.24) 



where we have introduced the functions 



f±(t) = ^{cospt/i^ax, . . . , a n )} cos[2th 2 (a u . . . , a n )] <g) 2^ V / B }, (2.25) 



and 



/,(*) = ^{cospt/i^ai, • • • , a„)] ® 2-^/b}, 



(2.26) 



where hi(ai, a 2 , • • • , <x n ) and Wai, «2, • • • , a„) are 



/ii(«i,a 2 , ...,«„) 



h 2 {a 1 ,a 2 ,. ■ -,a n ) 
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J2 a l J + J -> 
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\ 



E a 5 J - J +' 



Calculating the traces over the spin bath in the eigenbasis of J3 and J 2 using 



We find 



and 



J±J*\j,rn,x) = (j±m)(j^m+ l)\j,m,x). 



M)4iwi ,a( *g^ 



i=l ji,mj 



/.(«) = (f[ETW,A)]^f, 



i=l ji.mt 



here,C and 77 denoted are as 



and 



and also, we have 



c = £#«?, 



i=l 



rj^^hi a 2 , 
i=i 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



(2.31) 



where we have introduced the quantity h(j, ±m) = (j ± m)(j ^f m + 1). 

Thus we have determined the exact dynamics of the reduced system: The density matrix 
Ps(t) of the central spin is given through the components of the Bloch vector which are pro- 
vided by the relations (23), (24) and (30), (31). We note that the dynamics can be expressed 
completely through only two real- valued function f±(t) and f z (t). This fact is connected to 
the rotational symmetry of the system. Also from the overall role of coupling coefficients and 
degeneracy coefficients in the relations (30) and (31), it will be shown that the study of control 
over decoherence is on coupling coefficients and degeneracy coefficients that you'll see Sec. (3). 
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3 Example 

In this example, we consider the central spin with two different bath by two different coupling 
coefficients. This means we consider two layers with different radii. The choice of different 
coupling coefficients certainly will affect on degeneracy coefficient because the behavior of 
particles in each layer is different from other layers and this difference is to express the con- 
sidering by different coupling coefficients and different degeneracy coefficients. So,we consider 
our hamiltonian as 

H = 2(a + E_ + a_E + ), (3.32) 

where we define H as 



'+■ 



2 

E + = E a » J + 
11=1 

2 

According to the introduced model, we can express two real- valued functions f±(t) and f z (t) 



as 



f±{t) = 2^ 1^ ^{N 1 ,j 1 )T{N 2 ,j2) ^TJv^ — ' ( 3 - 33 ) 



Jl,mi J2,rn,2 

and 

m = E E n^nm^nf-^^l, (3 .34) 

where j3 and 7 are denoted as 

/3 = a\h(ji, mi) + alh(j 2 , m 2 ), 

and 

7 = alh(jt, -mi) + a%h(j 2 , -m 2 ). 

The explicit solution constructed in the previous section takes on a relatively simple form in the 
limit N — y 00 of an infinite number of bath spins [10] . Because N is Ni + N 2 , we should discuss 
about Ni — > 00 or N 2 — y 00 or both, in Ref.JTD] for large N, the corresponding correlation 
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function is obtained. But here by considering two layers, we obtain the states of N± — > oo and 
AT 2 — y oo for corresponding correlation functions as follows: 

\ 
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(3.35) 



(3.36) 



Of course, we assume a non-trivial finite limit iVj — > oo, therefore we rescale the coupling 
constant as [TO], 



O; 



a,- -)■ 



/]V7 



(3.37) 



Using this approximation in Eq.(19), we can rewrite f± and f z as 

\ 
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(3.38) 



(3.39) 



Z=0 fc=0 

The up state is a state with consideration to N\ — > oo and A^ — >• oo, which is for the most 
general level cases non-trivial for an infinite limit. But there might be a layer which has a 
finite number of particle and another one with an infinite number of particle (of course here, 
it is understandable that infinite means the order of Avogadro's number because this size of 
particles in the scale of our study, accounts for infinity). We assume that N 2 is limited and Ni 
is unlimited. Of course in computing because of the symmetry of the system, there isn't any 
difference is taking A^ limited and N 2 unlimited with the previous case. So, for Ai — > oo we 
can obtain as 

_. oo oo 

/±(*) = T^Yl S S T ( iV2 '^ 2 ) cos(2ta 2 ^h(j 2 , m 2 )) cos(2ta 2 ^/h(j 2} -m 2 )) 



J2,iri2 n=0 k=0 



(2w)! ' (2fc)! 

and 
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x ^^!^!^«? n+fc) > (3-40) 



^(*) = i S T ( JV 2, J2) cos(4ta 2 V^(j 2) m 2 )){l - 2 v / 2to 1 DF[ v / 2to 1 ]}. (3.41) 

Note that Dawson Function is closely related to the error function erf, as 



DF(x) = exp(— x )erfi(x), 

where erfi is the imaginary error function, erfi(x) = —ierf(ix) 



4 approximation techniques 

In this section we will apply different approximation techniques to the spin star model intro- 
duced and discussed in the previous section. Due to the simplicity of this model we can not 
only integrate exactly the reduced system dynamics, but also construct explicitly the various 
master equations for the density matrix of the central spin and analyze and compare their per- 
turbation expansions. In the following discussion we will stick to the Bloch vector notation. 
Each of the master equations obtained can easily be transformed into an equation involving 
Lindblad superoperators using the translations rules 

1 + UJ ^3 , , (A AO\ 

Ps = ^ 1- u+<J- + UJ-<J+. (4-42) 

The second order approximation of the master equation for the reduced system is usually ob- 
tained within the Born approximation [5]. It is equivalent to the second order of the Nakajima- 
Zwanzing projection operator technique. In our model the Born approximation leads to the 
master equation 

Ps(t) = ~ dstr B {[H, [H, p s (s) <g> Pb(0)]]} = -8fti / ds(u(s) 3 a 3 + cu(s) + a_ + w(s)_o+), 
Jo Jo 

(4.43) 
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where the bath correlation function is found to be 

fii = ^tr B {E + E_} = —tr B {( ai Jl + a 2 J 2 + )( ai J l _ + a 2 J 2 )} = -(c£/Vi + a 2 2 iV 2 ). (4.44) 

It is important to notice that Qi, as well as all other bath correlation functions are independent 
of time. This is to be contrasted to those situations in which the bath correlation function 
decay rapidly and which allow the derivation of a Markovian master equation. The time- 
independence of the correlation function is the main reason for the non-Markovian behavior 
of the spin bath model. The integro-differential Eq.(43) can easily be solved by a Laplace 
transformation with the solution 

f ± (t) = ^- = cos(2tV~S), (4.45) 

fM = ^cos(2tV25), (4.46) 

where 6 is denoted as 

5 = a 2 JVx + a\N 2 . (4.47) 

In many physical applications the integration of the integro-differential equation is much more 
complicated and one tries to approximate the dynamics through a master equation which is 
local in time. To this end, the terms oo±(t) and ui 3 (t) under the integral in Eq.(43) are replaced 
by u±(t) and uj 3 (t), respectively. We thus arrive at the time-local master equation 

d 



-p a (t) = -45 / ds(u 3 (t)a 3 + u+(t)a- + u-(t)<r+) = -A5t(u 3 (t)a 3 + cu+(t)<j- + w-(t)a+), 
at Jo 

(4.48) 



which is sometimes referred to as Redfield equation. Also this master equation is easily 
solved to give the expressions 

f ± (t) = exp(-2St), (4.49) 

f z (t)=exp(-A5t). (4.50) 



14 



1.0 



0.5 



n+i 0.0 



-0.5 
-1.0 



\ '^"^ 










-- a 1= &3 











0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 2: The equal spin numbers for correlation function f± [see Eqs.(33)], with a = 0.1 and 
Ni = 20. 
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Figure 3: The changes of spin number of layers for correlation function f± [see Eqs.(33)], with 
a = 0.1, JVi = 20 and N 2 = 100. 



The Redfield equation is equivalent to the second order of the time-convolutionless projection 
operator technique, finally, In order to obtain, a Markovian master equation, i.e. a time-local 
equation involving a time independent generator, one pushes the upper limit of the integral 
in Eq.(48) to infinity, as other studies of the master equation. This limit leads to the Born- 
Markov approximation of the reduced dynamics. In the present model, however, it is not 
possible to perform this approximation because the integrand does not vanish for large t. 
Thus,the Born-Markov limit does not exist for the spin bath model investigated here and the 
description of decoherence processes requires the usage of non-Markovian methods. 
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Figure 4: The equal spin numbers for correlation function fs [see Eqs.(34)],with a = 0.1 and 
N 4 = 20 . 
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Figure 5: The changes of spin number of layers for correlation function f^ [see Eqs.(34)], with 
a = 0.1, Ni = 20 and N 2 = 100. 
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5 conclusion 



With the help of a simple analytically solvable model of a spin star system, we have considered 
layered environment model with one-qubit in centr of layers and every layer is constructed by 
some spins. Of course, we have assumed these layers have different radii which means different 
coupling coefficients. By considering the Fig. (2) and Fig. (4), and selecting equal degeneracy 
coefficients of both layers and considering the influence of coupling coefficients in decoherence- 
time control it can be deduced that the fluctuations in f z is more intense than f±. Meanwhile 
variations of coupling coefficients, makes the fluctuations of f z greater than f±. In Fig. (4) 
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Figure 6: Comparison of limit iV — > oo with limit Ni — > oo or N 2 — > oo for the changes of 
correlation function f± [see Eqs.(38) and Eqs.(40)], with a = 0.1 and n = 20. 



one can see that the increase of fluctuations can allow us to control, the decoherence. but 
it's not good at f± because near zero, fluctuations are gentle (this means that the amplitude 
of oscillation is shorter and the course is slowly changing). Previous time of the system (at 
short period of times) is longer than the decoherence time at f z and is pulsed and is shorter 
at f± and continuous. As a result, in general case for f z and f± to increasing the power of 
coupling coefficient, yields greater control on decoherence and increase of decoherence time. 
Also, with by considering Fig. (3) and Fig. (5) and selecting inequal degeneracy coefficients, we 
find that the same effects of selecting equal degeneracy coefficients in Fig. (2) and Fig. (4) take 
place, but in this part, fluctuations are more and the carve shifts towards the vertical axis. It 
should again be emphasized that at f z transformation from coherence mode to decoherence 
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Figure 7: Comparison of limit iV — > 00 with limit Aq — > 00 or N 2 — > 00 for the changes of 
correlation function f 3 [see Eqs.(39) and Eqs.(41)], with a = 0.1 and n = 20. 



mode and vice versa is pulsed and at f± the decoherence is controlled continuously. Then, we 
have considered the limit iV — > 00 and for a special example (two layers with different couple 
coefficients) showed that N% — > 00 with N 2 — > 00 have equal result. By considering Fig. (6), the 
first chart is related to when Ni — > 00 and A^ — > 00, and the second chart is when Aq — > 00 
or ^2 —?• 00. In the first case we can find that the fluctuations are increased by the increase 
of coupling strength, this can also be seen in the second case too. But the difference between 
first and second situation is in decoherence-time control that in the first situation, this work is 
well done because the curve in the first graph is distant from the limit of zero and this means 
the increase in decoherence-time control. Thus the best case is when Aq — > 00 and N 2 — > 00 
occur simultaneously, which is more general and real. It is much better because decoherence 
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time is controled better. In Fig. (7) we have the same conclusion. In Fig. (6) and Fig. (7) we can 
see that we have more intensity in fluctuations of f z than f±, be side decoherence time at f± 
is better than f z because at f% transformation of states from coherence mode to decoherence 
mode and vice versa is pulsed. 

References 

[1] M.A.Nielsen and I.L.Chuang, Quantum Computation and Quantum Information (Cam- 
bridge University Press, Cambridge, 2000). 

[2] D.Loss and D.P.DiVincenzo, Quantum computation with quantum dots , Phys. Rev. A57, 
120 (1998). 

[3] B.E.Kane, A silicon-based nuclear spin quantum computer , Nature 393, 133 (1998). 

[4] D.D.Awschalom, N.Samarth, and D.Loss, Semiconductor Spintronics and Quantum Com- 
putation (Springer, New York, 2002). 

[5] H.P.Beruer, F.Petruccione, The Theory Of Open Quantum System (Oxford University 
Press, Oxford, 2002). 

[6] E.C.G.Sudarshan, P.M.Mathews and J.Rau, Stochastic dynamics of quantum-mechanical 
systems , Phys. Rev. 121, 920 (1961). 

[7] A.Hutton, S.Bose, Mediated entanglement and correlations in a star network of interacting 
spins , Phys. Rev. A 69 042312 (2004). 

[8] Y.Hamdouni, M.Fannes and F.Petruccione, Exact dynamics of a two-qubit system in a 
spin star environment, Phys. Rev. B 73, 245323 (2006). 

[9] Y.Hamdouni, An exactly solvable model for the dynamics of two spin-| particles embedded 
in separate spin star environments, J. Phys. A: Math. Theor. 42 315301 (2009). 



19 

[10] H.P.Breuer, D.Burgerth, F.Petruccione, Non-Markovian dynamics in a spin star system: 
Exact solution and approximation techniques , Phys. Rev. B 70 045323 (2004). 

[11] V.Semin, I.Sinayskiy and F.Petruccione, Initial correlation in a system of a spin coupled 
to a spin bath through an intermediate spin , Phys. Rev. A. 86, 062114 (2012). 

[12] H.Krovi, O.Oreshkov, M.Ryazanov, D.A.Lidar, Non-Markovian dynamics of a qubit cou- 
pled to an Ising spin bath , Phys. Rev. A 76(5), 052117 (2007). 

[13] E.Ferraro, H.P.Breuer, A.Napoli, M.A.Jivulescu, A.Messina, Non-Markovian dynamics of 
a single electron spin coupled to a nuclear spin bath , Phys. Rev. B 78(6), 064309 (2008). 

[14] D.Rossini, P.Facchi, R.Fazio, G.Florio, D.A.Lidar, S.Pascazio, F.Plastina, P.Zanardi, 
Bang-Bang control of a qubit coupled to a quantum critical spin bath , Phys. Rev. A 
77(5), 052112 (2008) 

[15] C.M.Dawson, H.P.Hines, R.H.Mckenzie, G.J.Milbrun, Entanglement sharing and deco- 
herence in the spin-bath , Phys. Rev. A 71 052321 (2005). 

[16] W.H.Zurek. Decoherence and the transition from quantum to classical REVISITED. 



arXiv:quant-ph/0306072, June 2003 



[17] E.H.Lieb, T.D.Shultz, D.C.Mattis, Ann. Phys. (N.Y.) 16 407 (1961) also reprinted in The 
Many-Body Problem by D. C. Mattis (World Scientific, Singapore, 1993). 

[18] A.Imamoglu, D.D.Awschalom, G.Burkard, D.P.Divincenzo, D.Loss, M.Sherwin, A. Small, 
Quantum Information Processing Using Quantum Dot Spins and Cavity QED , Phys. Rev. 
Lett. 83 4204 (1999). 

[19] A.Hutton, S.Bose, Mediated Entanglement And Correlations In A Star Network Of In- 



teracting Spins , quant-ph/ 0208114 (2002). 



20 

[20] J.Wesenberg, K.M0lmer, Mixed collective states of many spins , Phys. Rev. A 65 062304 
(2002). 



